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Abstract 

The low temperature heat capacity Cp of Fei-ajGa^^ ahoys with large magnetostriction has been 
investigated. The data were analyzed in the standard way using electron (7T) and phonon (/3T^) 
contributions. The Debye temperature decreases approximately linearly with increasing Ga 
concentration, consistent with previous resonant ultrasound measurements and measured phonon 
dispersion curves. Calculations Qo from lattice dynamical models and from measured elastic 
constants Cn, C12 and C44 are in agreement with the measured data. The linear coefficient of 
electronic specific heat 7 remains relatively constant as the Ga concentration increases, despite the 
fact that the magnetoelastic coupling increases. Band structure calculations show that this is due 
to the compensation of majority and minority spin states at the Fermi level. 

PACS numbers: 75.80.+q, 65.40.Ba, 62.20.Dc, 71.15.Mb 
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I. INTRODUCTION 



Fei-xGax alloys are known for their large magnetostriction. Values of the tetragonal mag- 
netostriction coefficient, |Aioo, can reach values as high as ~ 400 ppm for certain alloy com- 
positions and heat treatments.-'^'^ The sharp rise in the magnetostriction near 19at.%Ga 
composition can be explained by a simultaneously increasing magnetoelastic coupling, bi, 
and decreasing tetragonal shear modulus, C'.^ The decrease in C has been independently 
measured via resonant ultrasound techniques^'^'^ and neutron scattering''. However, the 
nature and characterization of the large increase in magnetoelastic coupling with composi- 
tion has been difficult to determine. It has been suggested that the increase in bi is related 
to short-range ordered clustering of the Ga atoms prior to the formation of long-range 
ordered structures near 19at.%Ga.- Below ~19%Ga, Fei_xGa.x alloys are disordered and 
crystallize in a body-centered-cubic (bcc) a-Fe (A2) structure. Above this composition, two 
ordered phases are possible; DO3 and B2. For itinerant magnetic alloys without significant 
short-range ordering (local) effects, bi depends on the spin-orbit coupling of electrons near 
the Fermi level. In strongly ferromagnetic alloys, the compositional dependence of the mi- 
nority spin electronic density-of-states (DOS) at the Fermi level, nilsp), can be related to 
magnetostriction.- Measurements of the linear coefficient of the electronic specific heat at 
low temperatures, 7, (also called the Sommerfeld constant) are directly proportional to the 
total electronic DOS at the Fermi level, n{eF), and can be used to characterize the ori- 
gins of magnetostriction. This is most clearly demonstrated in Ni-Fe alloys, where the zero 
magnetostriction composition (permalloy) corresponds to a full majority spin band and a 
minimum in nilep) leading to a minimum in 7, as predicted in the "split-band model" of 
Berger et al. (Ref.-). We have undertaken a study of the low tempeature heat capacity 
of Fe-Ga alloys to determine if variations in 7 are present that can be correlated with the 
large increase in magnetostriction. In addition, the lattice contribution to the specific heat, 
characterized by the Debye temperature 6d, indicates strong lattice softening with added 
Ga in agreement with ultrasound and neutron scattering data. 
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II. SAMPLE PREPARATION 



Single crystal alloys of Fei^^Gax were grown by the Bridgman technique (see Ref.- for 
more details of sample preparation). Gallium (99.999% pure) and electrolytic iron (99.99% 
pure) were cleaned and arc melted together several times under an argon atmosphere. To 
prepare single crystal samples, the as-cast ingot was placed in an alumina crucible and heated 
under a vacuum to 1500 C. After reaching 1500 C, the growth chamber was backfilled with 
ultra high purity argon to a pressure of 2.76 x 10^ Pa. Following pressurization, heating 
was continued until the ingot reached a temperature of 1600 C and held for 1 hour before 
being withdrawn from the furnace at a rate of 5mm/hr. Following crystal growth, the 
ingot was annealed at 1000 C for 168 hours (using heating and cooling rates of 10 degrees 
per minute). Small parallelepipeds (2mm xlmm x 0.5mm) were cut from the ingot by 
wire electrical discharge machining and cleaned by acid etching and polished on one side. 
Samples were sealed in a quartz tube and annealed at 1000 C for 4 hours and furnace cooled 
down to room temperature. Composition measurements were done by energy-dispersive 
spectrometers (EDS) in a JEOL 840A Scanning Electron Microscope (SEM). 

III. MEASUREMENTS 

Heat capacity measurements were carried out using a Quantum Design Physical Prop- 
erty Measurement System. The addenda were measured separately immediately before the 
sample measurement and subsequently subtracted. Heat capacity Cp data for slow-cooled 
Fei_i.Gaa; alloys are shown in Fig. [H 

IV. DATA ANALYSIS 

At low T, the heat capacity for a soft ferromagnet is given by the formula 

Cp = jT + pT^ + aTl, (1) 

where the terms represent the electronic, phonon, and spin-wave contributions, respectively. 
The spin wave contribution has been measured at very low temperatures and is estimated 
to be about a factor of fifty times smaller than 13.—^^ Our own analysis of the pure a-Fe 
{x = 0) data demonstrates that the heat capacity is not sensitive to the small spin wave 
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FIG. 1: Heat capacity Cp versus temperature T data for slow-cooled Fei-^Gaa; alloys. 



TABLE I: Parameters obtained from fits by Eq. (2) to slow-cooled Cp/T data. 

^t-%Ga7(;^)/3(;^) Bz. (K) 
0.0 4.89(1) 0.0187(1) 470.2(8) 
9.0 5.03(1) 0.0276(2) 413(3) 

15.4 4.46(1) 0.0332(1) 388(3) 

17.5 4.69(1) 0.0582(3) 322 
19.4 5.21(1) 0.0521(2) 334(4) 

term proportional to a below 10 K. In this limit, the spin wave term is ignored and the heat 
capacity can be written in the following form 

^ = 7 + /5Tl (2) 

The plot of the data as Cp/T vs. is shown in Fig. ([2]). The Sommerfeld constant 
(y- intercept) and lattice specific heat coefficient (slope) can then be obtained by a linear 
least-squares fit to the plot of Cp/T vs. T^. The parameters obtained from the fits are shown 
in Table [H The Debye temperature, Qd can be derived from (3 in the procedure described 
below and is also shown in Table [H The value obtained for 7 and G d for pure a-Fe are 
consistent with literature values.— »^»^»^»^»^ 
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FIG. 2: Heat capacity divided by temperature Cp/T versus for sfow-cooled Fci-x^Ga^ aUoys. 



A. Electronic heat capacity 

Figure [3] shows that the electronic coefficient 7 remains essentially constant as the Ga 
concentration increases. At low temperatures, the electronic heat capacity is proportional 



the universal gas constant. 

To better understand our experimental results, we also performed density functional 
calculations for Yei_xG&x alloys, using the highly precise full potential linearized augment 
plane wave (FLAPW) method. No shape approximation is assumed for the charge, poten- 
tial, and wave function expansions in the entire space. We used the generalized gradient 
approximation^^ for the description of the exchange correlation interaction. The conver- 
gence against parameters such as the number of k-points and energy cutoff was carefully 
monitored. We used a (2 x 2 x 2) supercell that comprises 16 atoms throughout the cal- 
culations and we studied cases with x = 0.0 (pure bcc a— Fe), x = 0.0625 (FeisGai), x = 
0.125 (Fei4Ga2) and x = 0.1875 (FeisGas). For the Fei4Ga2 and FeisGas cells, there are 
several different ways to arrange Ga atoms on the bcc lattice sites are related to different 
short-range ordered structures. Results reported below correspond to their minimum energy 
configurations. Fig. Hlshows the calculated electronic DOS for Fei^xGa^ alloys with x = 0.0, 
0.0625, 0.125 and 0.1875, respectively. As x increases, n^ep) changes oppositely in the two 
spin channels. In the majority spin part, the Fe ci-holes are gradually purged and n| (ep) 



to the DOS at the Fermi level n{eF), according to the formula 7 = 



n{eF), where R is 
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FIG. 3: (a) The electronic coefficient of the heat capacity (7) of slow-cooled Fei-xGa^ alloys 
plotted as a function of Ga concentration, x; experimental data (empty circles) and calculated 
(filled circles). In experimental data, the statistical error bars are smaller than the symbol size, 
(b) The spin-projected electronic density of states of majority (up triangles, solid line) and minority 
(down triangles, dotted line) at the Fermi level n^ep). 



falls monotonically. For x =0.15, the Fe-d band in the majority spin channel is completely 
filled. Meanwhile, the number of non-bonding states around the Fermi level grows steadily in 
the minority spin channel. The trend of the n (ep) versus x curve reasonably matches with 
the experimental data of 7, as shown in Figl31 In the calculation, the small dip around x = 
0.15 mainly stems from the elimination of the Fe majority spin rf-holes due to the presence 
of Ga atoms. This dip feature is seen also in the experimental data, however it is a small 
effect that is on the limit of the sensitivity of the technique. 
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FIG. 4: The calculated electronic density-of-states for Fei_2,.Gax. alloys with black solid lines for 
X = (i.e., pure bcc Fe), green dashed-dotted lines for x = 0.0625, red dashed lines for x = 0.125, 
and blue dotted line for x = 0.1875. The positive side is for the majority spin channel while the 
negative side is for the minority spin channel. Zero energy is for the position of the Fermi level. 

B. Lattice heat capacity 

For a single, isotropic phonon mode with sound velocity c, the contribution to the low- 
temperature molar heat capacity in the Debye model is 

2n^RklV , 

- -^w^^ ' 

where V is the volume of the primitive cell, R is the gas constant, and h is Planck's con- 
stant divided by 2tt. For a general cubic crystal with elastic anisotropy and three phonon 
polarizations, the form of the heat capacity is the same as Eq. ([3]), with c replaced by the 
effective sound velocity c, which is obtained by averaging the inverse-cubed sound velocities 
over all possible propagation directions and modes 




(4) 
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FIG. 5: The Debye temperature Qd plotted as a function of Ga concentration x as determined 
from the data are shown as sohd points (•). The solid line is a linear fit to the slow-cooled data. 
Various calculated 0d values are indicated by the open circles (o), as referenced in the text. 



where i>i{6,ip) is the sound velocity in a crystalline direction given by 6,f. The Debye 
temperature is defined as 

e.^l^m*. (5) 



kB\V ^ 

such that the low-temperature heat capacity can be written 

Cv - ^ (6) 

and the slope of CyjT versus T2 equals /? = (1943.9 J mol-^K-i)e^l Experimentally de- 
termined values of are given in Table [T] and are also shown in Fig. [51 G/j decreases 
approximately linearly with increasing Ga concentration, consistent with previous measure- 
ments of the tetragonal shear modulus C via resonant ultrasound^ and neutron scattering^. 

In an effort to affirm the validity of the measured Debye temperatures, we numerically 
calculated via three methods outlined below. These various estimates of G^) are also 
shown in fig. [5l 

(1). Integral: The simplest way to calculate the Debye temperature from a lattice 
dynamical model is to calculate the heat capacity by integration. Given the model lattice 
dynamical parameters (obtained from fits to inelastic neutron scattering data, for example), 
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FIG. 6: (a) Phonon density-of-states g{E) versus energy E for a-Fe as determined from a force 
constant model, (b) g{E)/E'^ versus E showing the hmiting value (dashed line) as E approaches 
OmeV. (c) Cy/T versus T for the data (empty circles), and as calculated from Eqn. ([7]) (solid 
line). The best linear fit to the calculation (dashed line) gives a Debye temperature of 464 K. 



the phonon density-of-states, g{E), can be calculated. The phonon DOS can then be used 
to calculate the heat capacity according to 

2 



Cy = 3R 



9iE) 



E 



(7) 



e*sT - 1 

At low temperatures, the result of this calculation can be fit to the limiting form in Eq. ([6]) 



to obtain Qd in the same way as it was obtained from the measured data above. Fig. EJ^a) 
shows the calculated DOS for a-Fe where the force constants are obtained by fitting the room 
temperature phonon dispersion curves measured by inelastic neutron scattering.— Using this 
DOS, we can calculate the heat capacity at low temperatures via Eq. ([7]) and, subsequently, 
fit the curve up to 10 K as shown in Fig. [Hl^c). The fitted slope, 1.94x 10^^ Jmol^^K"^, can 
then be used to obtain Od = 464 K via Eq. (Q. 

(2). Extrapolation: One can obtain Qd directly from the DOS without calculating the 
heat capacity. The dependence arises from the low energy quadratic energy dependence 
of the DOS (resulting from a linear dispersion relation). Plotting YmvE^o g{E) / E"^ gives the 
coefficient of the low energy quadratic form of the DOS, which leads directly to the average 
c: 



^ ^.m 



(8) 



27r2 \e^o E2 

Then, Eq. ^ can be used to calculate 9^. Fig. [6](b) shows g{E)/E'^ as calculated for 
a-Fe. The (average) zero frequency limit of this function has a value of approximately 
4.7x10"^ meV"^ Using Eqs. ([5]) and ([8]), Od is calculated to be 463 K. 

(3) Elastic constants: Additionally, the Debye temperature can be calculated directly 
from the elastic constants. Any cubic system can be described by three independent elastic 
constants; Cu, Cu and C44. Generally these are experimentally measured as two transverse 
modes, C44 (rhombohedral shear) and C = l/2(Cii — Cu) (tetragonal shear), and one 
longitudinal mode, K = l/3(Cii + 2C12) (bulk modulus). Given these elastic constants, 
the sound velocities z/i(6', 0) can be calculated for any propagation direction {6, 0) using the 
Green-Cristoffel equations: 



(Cu - C44)ql + Cmq'^ - puj'^ (Cu + Cii)q^qy (C12 + Ci4)q^q., 

(C12 + Cii)qyqa, [Cii - Cii)ql + d^q^ - pu"^ {Cu + C44)qyqz 

{Cu + C44)qzqx {Cu + ^44 

(9) 

Solutions of Eqs. (Q yield the sound velocities ^'^(q) = ujj{q)/q for each phonon branch 
along the crystal direction given by wave vector q. The average sound velocity can then 
be evaluated for an anisotropic cubic crystal by numerical averaging as shown in Eq. (jlj). 
Subsequently, Eq. ([5]) can be used to calculate Qd- 
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For a-Fe, the zero-Kelvin elastic constants are extrapolated to be Cu = 243.1 GPa, Cu = 
138.1 GPa and C44 = 121.9 GPa.— Calculation of the Debye temperature based on the above 
method gives Qd = 478 K which compares favorably to the value determined from fits to 
the low temperature calorimetry data, as shown in Fig. O This can also be compared to the 
result obtained via the de Launay formula; a general semi-analytic function to determine 
Qd{T) for cubic metals from the elastic constants.— i^^i^i Given the elastic constants listed 
above, this formula yields Qd{T = OK) = 477K.— 

The elastic constants for quenched Fegi.sGaig.y were independently measured by Clark 
et al. (Ref.^) as a function of temperature. These were extrapolated to OK and substi- 
tuted into the de Launay formula yielding 6£i(T = OK) = 386 K. They were also used via 
Eqs. dnHHandE]) to calculate Qr, = 396 K. As can be seen in Fig.O these values are generally 
consistent with trends in the data. 

V. DISCUSSION 

We find only a weak correlation of 7 and Aioo with composition. Band structure calcu- 
lations within density functional theory show that the small change in 7 is caused by the 
cancellation of two effects; a simultaneous depletion of holes in the majority spin band and 
the expected increase in riilep) in the minority band. Thus, while the heat capacity is not a 
sensitive probe of the relevant electronic states for magnetostriction in the case of Fei_a;Gaa;, 
the agreement with band structure calculations is a positive step. The increase in n{eF) 
after x > 0.15 manifests the presence of non-bonding states around the Fermi level in the 
minority spin channel. As revealed in our previous studies,- these states play a key role for 
the increase of Aqoi and the role of the effect of short-range ordering is an open question. 
Of course, more analysis on the details of wave functions, such as their magnetic quantum 
numbers, are needed for correct prediction. Results concerning the theoretical prediction of 
Aooi will be published elsewhere. 

In addition to the electronic behavior, the low temperature heat capacity can indepen- 
dently probe the lattice softening by determination of the Debye temperature. We find that 
the decrease of the Debye temperature agrees with previous estimates of lattice softening. 
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VI. SUMMARY 



We have measured the heat capacity as a function of temperature for Fei-ajGa^, 0.0 
< X < 0.194 sohd solutions crystalhzing in the bcc structure. The Debye temperatures 
follow a linearly decreasing trend with increasing Ga concentration consistent with known 
lattice softening. The electronic coefficient of the specific heat remains relatively constant, 
in agreement with band structure calculations. 
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